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Abstract 

We calculate the yields of a variety of hadrons for RHIC and LHC energies 
assuming thermodynamical equilibration of the produced minijets, and using as 
input results from pQCD for the energy densities at midrapidity. In the calculation 
of the production of partons and of transverse energy one has to account for nuclear 
shadowing. By using two parametrizations for the gluon shadowing one derives 
energy densities differing strongly in magnitude. In this publication we link those 
perturbatively calculated energy densities of partons via entropy conservation in 
an ideal fluid to the hadron multiplicities at chemical freeze-out. 



1 



1 Introduction 



Particle production in high-energy heavy-ion reactions at the BNL-RHIC 
and CERN-LHC colliders will soon provide interesting insight into nuclear 
modifications of semi-hard processes 0. This is because pQCD processes in- 
volving gluons in the initial state will dominate the inelastic AA cross-section 
at collider energies. In particular, they might lead to a better understanding 
of the gluon distribution in large nuclei, which is not accessible in DIS. 

In [0] the effect of nuclear shadowing of the parton distribution functions 
on the charged particle multiplicity at midrapidity has been investigated as- 
suming no rescattering between the produced minijets and the hadrons they 
fragment into. Here, we will take the opposite point of view and assume max- 
imal rescattering, i.e. local thermal and chemical equilibrium of the minijets. 
We compute final-state hadron multiplicities of various hadron species under 
the assumption of entropy conservation. 

In |J we calculated the initial conditions at RHIC and LHC by means 
of pQCD above the semihard scale px — 2 GeV to derive the number and 
energy densities of partons at midrapidity. In that calculation we explicitely 
included the shadowing effect on the parton distribution functions entering 
the formulas for the production of flavor / = g, q, q in the minijet approach. 
We employed two different parametrizations for the shadowing effect account- 
ing for weak and strong gluon shadowing, respectively, shown in figure |]. A 
direct consequence of the shadowing effect is the decrease in the production 
of partons of given momentum px, i.e. a decrease of transverse energy pro- 
duction at midrapidity (for a more detailed description of the calculation, 
the variables, and the hard partonic subprocesses see |3], f|, |5|, 0). 
We calculate the first moment with and without shadowed pdf's and 
with a cut-off function e(y) ensuring that we only count scatterings into the 
central rapidity region (\y\ < 0.5): 
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Figure 1: Shadowed parton distributions as parametrized in and in j^J. 





no shad. 


weak shad. 


strong shad. 


o-9 (E T ) 


26.74 


27.25 


17.80 


a q (E T ) 


3.20 


2.87 


1.00 


a~ q (E T ) 


1.89 


1.73 


0.59 



Table 1: First ET-Moment in mb GeV at midrapidity (\y\ < 0.5) for par- 
ton production at RHIC for the different shadowing parametrizations and no 
shadowing, respectively. 

The results for the production of transverse energy by the semihard partons 
at RHIC for unshadowed, weakly shadowed, and strongly shadowed gluons 
are shown in table [I] (for the results from the different subprocesses see 0). 

By dividing the transverse energy E T ' (6 = 0) = T AA (0)a f hard (E T ) hard of 
parton species / by the initial volume 

Vi = 7iR 2 A Ay/p , R A = A 1 ' 3 x 1.1 fm (2) 

we found the energy densities for the three different cases as shown in table 
|2|. One sees that the results for the energy density with the weak shadowing 
are even slightly increased due to the onset of Rg(x,Q 2 ). We also included 
the contributions of the soft processes as discussed in ||. Motivated by 
the factorization in QCD one can assume that the production of transverse 
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no shad. 


weak shad. 


strong shad. 




60.0 


61.2 


40.0 


e q 


7.2 


6.4 


2.3 


Eg 


4.3 


3.9 


1.3 



Table 2: Energy densities for unshadowed, weakly and strongly shadowed 
gluons in GeV/fm 3 for Au + Au at RHIC. 

energy in AA collisions can be split up into a hard and a soft contribution as 
Mb) = T AA (b) [aZ rd (E T )Z rd + < P (Et)T] ■ (3) 

With an energy independent value of a PP = 32 mb one derives || a pp (E T ) p s p = 
15 mb GeV. With 

TauAu — 29/mb one can derive the soft contribution 
(i.e. the one for p T < 2 GeV) to the energy density for RHIC as e so f t = 
33.7 GeV/fni . If one assumes the soft contribution to be independent of 
the shadowing phenomenon one can derive the change in relative weight of 
hard to soft processes due to the shadowing of the hard contributions. For 
the total first Et moment of gluons, quarks, and antiquarks one finds that 

OLrd + a lard + °lard) ( E t) hard ec l uals 3L8 mb GeV for no, 19.4 mbGeV for 
strong, and 31.9 mbGeV for weak shadowing. 

When comparing to the soft contributions one finds with owd (Er)h = 
(a 9 + o q + o-^) hard (ET) hard that the ratio of soft to hard contribution 

p _ O~soft (E T ) S ( ,x 

&hard {^T/h 

is R s h = 0.47 for no, R s h = 0.77 for strong, and R s h = 0.47 for weak shadow- 
ing, respectively. This implies that at RHIC the soft component could even 
dominate if it were energy independent and also unaffected by the shadowing 
effect. 

The same analysis was done for LHC and we found for the first Et 
moments in the three cases the results depicted in table [| For LHC we 
only calculated the contribution of the gluons that strongly dominate all 
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no shad. 


weak shad. 


strong shad. 


a 9 (Et) 


513.01 


286.87 


60.39 



Table 3: First Ex-Moment in mb GeV at midrapidity \y\ < 0.5 for gluon 
production at LHC for the different shadowing parametrizations and no shad- 
owing, respectively. 



partonic processes due to the large distribution function at small momen- 
tum fractions. The results for the energy densities for no, strong, and weak 
gluon shadowing, then are e„ = 1229.7 GeV/mi , e g = 144.8 GeV/fni , and 
Eg = 678.6 GeV/fm 3 . With Tp^pb = 32/mb one can derive the energy density 
from the soft part and has £ S oft = 35.8 GeV/fm 3 which is slightly larger than 
at RHIC due to the larger nuclear overlap function, i.e. the larger number of 
effective scatterings in the Glauber picture leading to the transverse energy 
production. For the various shadowing parametrizations one finds for the 
total first Et moment {a 9 hard + a q hard + a q hard ) (E T ) hard is 513 mbGeV for no, 
60 mbGeV for strong, and 287 mbGeV for weak shadowing. The relative 
weight 

between soft and hard contributions therefore changes and becomes R s h = 
0.029 for no, R s h = 0.25 for strong, and R s h = 0.052 for weak shadowing, 
respectively. 

Therefore the soft contributions gain much more weight in this naive picture 
due to the strong effect of shadowing on the small-x gluons (for further details 
of the calculation, the shadowing parametrizations, etc. see 0). 



2 Hadron Multiplicities at Chemical Freeze- 
Out 

Having calculated the energy densities at midrapidity for RHIC and LHC in 
pQCD we connect £j to the number of hadrons at midrapidity by assuming an 
ideal fluid that is characterized by entropy conservation from the quark-gluon 
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plasma to the hadron gas, dSi/dy = dSf/dy ||. We relate the energy density 
of the quark-gluon plasma to its entropy density via the bag model equation 
of state jrUI. We account for u, d, s quarks (with masses m u = m d = 0, 
m s = 150 MeV), the antiquarks, and gluons. 

The total produced entropy dSi/dy is obtained from the entropy density 
at midrapidity as dSi/dy = ViSi with the initial volume of the central region 
Vi = 7iR\ri, which is numerically Vi = 12.9 fm 3 for Au + Au and Vi = 
13.4 fm 3 for Pb + Pb at b = with r = 0.1 fm/c. Since we assumed an ideal 
fluid the total entropy is conserved throughout the expansion until freeze- 
out which is chosen here to happen at a temperature Tpo = 160 MeV. For 
simplicity we furthermore assume vanishing chemical potentials in the central 
rapidity region, i.e. that all conserved currents are identically zero. If this 
were not true one would have to multiply by factors exp(ui/T) (in Boltzmann 
approximation) . 

The entropy density of the hadronic fluid is calculated assuming an ideal 
gas composed of all hadrons up to a rest-mass of 2 GeV. Their respective 
occupation numbers are given by Fermi-Dirac or Bose-Einstein distribution 
functions, respectively. Thus, T FO and dSf/dy = dSi/dy determine the 
multiplicity of each hadron species uniquely |TT], [L2]. Feeding from post 
freeze-out decays of heavier resonances is also taken into account. 

3 Results 

With the model outlined above and the energy densities derived above we 
calculated the number of a variety of hadrons at midrapidity. We also include 
the multiplicities due to the soft contributions and quote the initial temper- 
atures for a QGP of three flavors. For LHC we derived the yields shown in 
table |4] and for RHIC the ones in table [|. For the latter one, one can clearly 
see that there is no change in the hadron yield for weak shadowing and the 
unshadowed case, respectively, due to the almost identical energy density 
serving as an input for the calculation. 
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LHC 


7T + + 7T 


K 





P 


A 




n 


T 


no shad. 


2680 


478 


32.1 


91 


58 


9.3 


1.4 


881 MeV 


weak shad. 


1720 


306 


20.6 


58.3 


37.2 


5.9 


0.9 


760 MeV 


strong shad. 


538 


95.9 


6.5 


18.3 


11.6 


1.9 


0.3 


516 MeV 


soft contrib. 


187 


33 


2 


6 


4 


0.7 


0.1 





Table 4: Hadron yields at freeze-out with initial conditions from pQCD for 
LHC. The initial temperatures are for a three flavor quark-gluon plasma with 
two massless quarks and m s = 150 MeV. Since the soft processes do not 
significantly contribute to the temperature we did not calculate Tj for the soft 
partons. 



RHIC 


7T + + 7T 


K 





P 


A 




Q 


T 


no shad. 


316 


56.3 


3.8 


10.7 


6.8 


1.1 


0.2 


433 MeV 


weak shad. 


316 


56.3 


3.8 


10.7 


6.8 


1.1 


0.2 


433 MeV 


strong shad. 


217 


38.7 


2.6 


7.4 


4.7 


0.8 


0.1 


383 MeV 


soft contrib. 


179 


32 


2 


6 


4 


0.6 


0.1 





Table 5: Hadron yields at freeze- out with initial conditions from pQCD for 
RHIC. 
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4 Conclusions and outlook 



We computed the rapidity densities of a variety of hadrons based on the 
assumption of entropy conservation of an ideal fluid. The entropy densi- 
ties were derived from the energy densities which in turn were caluclated by 
means of pQCD ||. We find that in the limit that the minijets equilibrate 
locally the effect of shadowing on the hadron yield is not as large as on the 
pure partonic degrees of freedom. This can be seen e.g. in the ratio of en- 
ergy densities between unshadowed and strongly shadowed gluons at LHC, 
which is about a factor of 9, while the ratio of the hadron yields is only a 
factor of 5. Since vanishing net baryon and strangeness densities were as- 
sumed, the relative depletion of shadowed to unshadowed gluon distribution 
is independent of the particle species. For RHIC, even without shadowing 
we only get about 300 pions since the perturbative calculation, entering via 
the energy densities, was performed with the cut-off px = 2 GeV. Therefore, 
the soft contribution constitutes a significant part of the transverse energy 
and therefore of the particle multiplicities fl2 [. E.g. in the UrQMD model 



the total pion yield at midrapidity is ~ 1100 [13] for RHIC. Therefore with 



an energy density of e s = 33.7 GeV/fm 3 , which was in extracted in || from 
CERN data, and assumed to be energy independent, we get 316+179 pions 
and therefore get a significantly smaller multiplicity compared to calcula- 
tions in the hadronic cascade model UrQMD. The key to this discrepancies 
could stem from the assumed energy independence of the soft contribution. 
For a high energy process one can write the interaction as a current-current 
interaction as 

a (p T ) ir r r (J-^j # 1 ® / 3 ® / 2 + e> (J^ (6) 

where the f n are non-perturbative IR dominated matrixelement of (lowest) 
twist n and H depicts the hard part. One should interpret the soft processes 
at least in part to stem from higher twist effects. For an unpolarized process 
the next-to-leading twist is r-4. In general one can relate the matrix elements 



of higher twist to the well known parton distribution functions of r-2. In JEI 
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the twist-4 correlation functions in a nucleus were parametrized as 

Tq t g = X 2 A 1/3 f i/A (x, Q\ A) = \ 2 A^R q , g (x, Q 2 , A)f i/p {x, Q 2 ) (7) 

where A 2 ~ 0.05 — O.lGeV 2 . This immediately implies that the soft processes 
are not independent of the energy since x ~ 1/ y/s and therefore should also 
increase from CERN-SPS to BNL-RHIC / CERN-LHC. Equation © in turn 
implies that the soft processes also are affected by shadowing; and since the 
DGLAP evolution between Q ~ 1 GeV and Q ~ 2 GeV is rather slow, one 
should expect no change in the relative weight between hard and soft pro- 
cesses. As a result one has two competing effects: as the cm. energy increases 
the soft production of transverse energy increases due to the decrease in the 
momentum fraction x and eq. (0). On the other hand the shadowing ratio 
Rg(x,Q 2 ) suppresses those contributions as x decreases. This topic requires 
a more detailed analysis elsewhere. 

The interesting feature at RHIC is the dependence of the hadron yields on the 
onset of gluon shadowing which differs substantially in the two approaches 
we compared. Since for the weaker gluon shadowing given in J7| the energy 
density for RHIC is more or less the same as in the unshadowed case we see 
the same results for the particle numbers in the two cases. However, for the 
strong gluon shadowing we see a depletion to about 70% of the unshadowed 
value. The effect of shadowing on the hadron multiplicities at central rapid- 
ity is found to be weaker than on the newly produced partons. 
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